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APPROXIMATE SOLUTIONS OF THE RUNGE-KUTTA EQUATIONS 

SUMMARY 


In this article we are concerned with the determination of the para- 
meters which enter the higher-order Runge-Kutta processes which are used for 
the numerical solution of differential equations. As the order of the process 
increases, the number and the complexity of the system of equations which 
determine the parameters in the Runge-Kutta formulation increases; also, the 
number of the evaluations of the function which is being integrated increases 
rapidly with the order. Since the latter is frequently the most time consuming 
computational aspect of the integration problem, approximate Runge-Kutta 
formulations which reduce the number of these evaluations without impairing 
the efficiency of the process are acceptable. Some of the techniques for obtain- 
ing these parameters — and the nature of the approximations engendered — are 
examined here for particular orders as these involve general aspects which 
persist and are as pertinent in the higher-order approximations. 


THE RUNGE-KUTTA METHOD 


We briefly indicate the Runge-Kutta process which leads to an algebraic 
system of equations whose solutions give the parameters involved. 

Consider the initial value problem 

y* = f(x, y) y(x 0 )=y 0 • U) 

We define the sequence of N quantities 
ki = hf(x 0 , y 0 ) 

k 2 = hf(x 0 + d 2 h, y 0 + /? 21 k t ) 
k 3 = hf(x 0 + a 3 h, y 0 + ^ 31 kj + /3 32 k 2 ) 
k 4 = hf(x„ + a^h, y 0 + /? 41 kt + ^ tt k 2 +p a k 3 ) 


( 2 ) 



k 5 - hf(xo + Qigh, yo + /Jgjki + $52k 2 + /J^kj + 


n - i 

k N = hf ( X , + V' y » + .? i % k j) 1 

and the linear form 

n 

y(x 0 + h)«y 0 + Yj K (3) 

i - 1 

such that the expansion of this in powers of h agrees with the exact solution of 
equation (i) to a prescribed number of terms. 

It is the usual practice in determining the parameters in equations(2) 
to expand these functions in Taylor series about the initial point (x 0 , yo) so that, 

to a given order, equation (3) is identical with the Taylor expansion 

1^2 j 3 jj 

y(x 0 + h) =y 0 + hf(x 0 , y 0 ) +— F(x 0 , y 0 ) +— r (x 0 , y 0 ) + 

(4) 

which is the true solution of equation ( 1) . Here the superscripts denote total 
derivatives with respect to x, which in terms of the partial derivatives in- 
creases in complexity as the order increases. On equating the respective 
coefficients of the partial derivative in these expansions to a given order in h, 
we are led to the system of equations which determine the parameters in 
equations (2). 
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Kopal II] , using the classical Taylor series expansion, derives — and 
gives a number of solutions for — the Runge-Kutta equations through the fourth- 
order. Shanks [2, 3] and Butcher [4-6] , again basing their work on Taylor 
series expansions , show how the Runge-Kutta formulations may be made for 
higher orders. The writer 1 , in a basically new approach to this classic Runga- 
Kutta process , obviates these tedious Taylor expansions and, using a method of 
quadratures, shows how to derive a matrix equation which gives the Runge- 
Kutta equations for any order. 

As the order of the process increases the number and complexity of Ihe 
system of equations which define the parameters in equations (2) also in- 
creases; and the minimum number of derivative evaluations (or', see equation 
(1) , evaluations of the function f(x, y) ) required for an exact solution of the 
Runge-Kutta equations for a given order increases rapidly. 2 Since the most 
time consuming computational aspect is the evaluation of the functions in 
equations (2) , approximate solutions which reduce the number of these deriv- 
ative determinations is desirable if it is not attained through an unacceptable 
reduction in accuracy. 

Shanks [2, 3] has pioneered in making these approximations and gives 
values for parameters in equations [2] for orders through eight. 3 However, 
since only the numerical values of these constructions are given, informative 
material on the techniques for obtaining these parameters, the nature of the 
approximations and errors, etc. , has naturally been suppressed. It is for 
this reason that the detailed analysis is given here for the fifth and sixth-order 
cases; and the nature of the analysis, as well as the errors entailed, which are 
shown here , will frequently demonstrate only simpler aspects of the arithmet- 
ically more complex higher-order cases. 


^osen, J. S. : The Runge-Kutta Equations by Quadrature Methods. To 
be published. 

2 The relation between the order and the number of function evaluations, 
N, is complex; Butcher [7] deals with this problem. Thus, both Shanks [2] 
and Butcher [4] have shown that the fifth order requires at least six evaluations 
for an exact solution. 

3 He also gives some exact solutions; e. g. , he gives the parameters 
in equations [2] with an adequate value of N to constitute a true solution of the 
Runge-Kutta equations in the classical sense. For the fifth order Luther and 
Konen [8] and Luther [9] give a number of these classical Runge-Kutta for- 
mulas while Cassity [10] considers the general solution for this order. 
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THE FIFTH-ORDER RUNGE-KUTTA EQUATIONS IN MATRIX FORM 


It will be convenient to express the Runge-Kutta equations in a matrix 
form. In addition to the advantage of compactness in this form, generalizations 
may more easily be extended to higher orders. 


Two basic sets of relations between the parameters in the Runge-Kutta 
process may be set down. 4 The first is 


i 0 
1 a 2 


i a 


®! 


(coj u) 2 • • • ajfc.) I 1 a j 


or 
N N 


0 \ 
m 

m 


m , 

“N / 


- (1 i/2 i/3 ... . i/(m + 1)) 

(5) 


where (m + i) is the order and N is the number of function evaluations shown 
in equations (2) . The first five equations in Table I are given by equation ( 5) 
when m = 4. 


The second set of relations between the parameters in equations ( 2) is 
given by 


a. . = P. , 

l + i i + i, 




• + 0 


i + i, 


i= 1, 2, ..... (N - i) 

( 6 ) 


We will now give another matrix equation which gives the remainder of 
the Runge-Kutta equations for the fifth order shown in Table I. 


4 See footnote 1 , page 3 . 
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TABLE I. THE RUNGE-KUTTA EQUATIONS OF THE FIFTH 
ORDER WITH FIVE EVALUATIONS 


Wj + 01 2 + Olj 

+ CJ4 

+ 01 5 

1 

a*) 

CJ 2 af2 + OljOg 

+ W4CV4 

+ Ol gO! 5 

1 

2 

(2*) 

u 2 0 i\ + 0130:3 

+ co 4 c*4 

+ oi 5 a| 

1. 

3 

(3*) 

oj 2 otl + 0130:3 

+ 0140:4 

+ w 6 aj 

1 

4 

(4*) 

O^a:* + 0130:3 

+ OI4G: 4 

+ 01505 

i 

5 

(5*) 

^3 032^2 

+ 01 4 ( 042^2 + 043^3) 

+ ^5^062^ 2 + 053 a 3 + 054°U) 

1 

6 

(6*) 

<^3 032^2 

+ Ol4(/3 42 0'2 +0430:3) 

+ <^5(052^2 + 053^ 3 + 054<*5) 

J. 

12 

(7*) 

a, 3 032°'2 

+ 014(^420'! + 043°^) 

+ 01 5 ( 052 0?2 +053^3 + 0540:^) 

20 

(8*) 

^3 032 ^2^3 

+ 0)4(/3420:2 + 043^*3 )** 3 

+ oi 6 (/3 52 q: 2 +0530:3 + 0540:4)0^5 

_1 

8 

(9*) 

t ^3 032O'2 a 3 

+ 014(^420:2 + 043 Q; 3) a '4 

+ aj S^052 Qf 2 + 053°'3 + 054 Q ’4) Qf 5 = 

15 

(10*) 

0J 3 032^2^3 

+ oi 4 (/3 420:3 + 0430:3)0:4 

+ Ol B (0jaa? 2 + £5303 + ^540:4)0:5 

J_ 

10 

(11*) 

LL, 3 03 2 2 Qf 2 

+ 014(^420:2 + 043 Q, 3) 2 

+ Ol5(05 2 Of 2 + 0530: 3 +054O:4) 2 

J_ 

20 

(12*) 


^4032043^2 

+ 01 5 [/3 32 0 530:2 + (0 42 o:2 + 043 a 3)054] 

1 

24 

(13*) 


w 4032 043 Q! 2^ a 3 + <*4) 

+ Ol 5 [/3 3 2 0530:2(^3 + “5) + (042^2 + 043°' 3) 054(^4 + O' 5)] - 

7 

120 

(14*) 


014/832 043O: 2 

+ Oi 5 [032053Of| +(0420:1 + 0430:1) 0 54 ] 

A 

60 

(15*) 



015032 043 054^2 

1 

120 

(16*) 


a 2 = 021 
a 2 = 031 + 032 
a 4 = 041 + 042 + 043 
a 5 ~ 051 + 052 + 053 + 054 
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Let us define the matrix 



o> 4 

CO 40 ; 4 


(o) 


(i) 

CO 4 C 4 

(°) 

74 

74 0) O!4 

rl l) 

75° ) / 3 54 



W 5 C 5 

7j 0) 

75 0) «5 

( 1 ) 

75 


( 7 ) 


where 


and 


/ n n n \ 

« 4 a 4 ... w N a N J 


032 0 0 . . 

£42 £43 0 

0 52 053 054 0 


B 8 B B 
P N2 P N3 P N4 P N5 


/ (n) (a) 

= 1,72 73 


0 


\ 


N 


■■ 

*N-i 0 / 


7^) , n = 0, 1 


(8) 


c s (,) 

~ 032 “2 




ci J) 

= 042 <*2 

+ 043 Q!3 


(9) 

;o) 

C N 

= 8 Qjj 

P fc2 2 

+ B 0<j + . 
P N3 3 

, . + 8 Oi Z 

P N,N-i N-l 



j = 1» 2, 


, m-i . 
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The matrix (7) then becomes 


0 = 


/ w 3 

“4 

^5 

/ ^ 3^3 

w 4 a 4 

« 5«5 

f ^*“3 

^ 4«4 


W 3 ^ 32«2 

a) 4 (^ 42 0! 2 + £430:3) 

a ’s(^ 52 Q! 2 + 

(a; 4^ 43 + CO5P53) 

(«g^ 54 ) 

0 

( Gl? 4/3 43 + 6^5^53) a 3 

(W 5 / 354 )Q !4 

0 

Y (<^ 4 ^ 43^4 + ^5^53^5) 

(Wg^ 54)0:5 

0 

\ (^ 5 ^ 54)^43 

0 

0 


( 10 ) 


Let us also define the matrix whose elements are given by ( 9) 


C = 


&32 a 2 P& a 2 P$2 a 2 

042 a 2 + 043 a 3 Z 3 42 a 2 + Z 3 ^! Pu a 2 + Z 3 #^ 

^52^2 + PsS&2 + ^54^4 Z 3 52 0! 2 + ^SS^I + 054 a 4 P® a 2 + ^53 a 3 + / 3 54 Q; I> 


Then (See footnote 1, page 3), 



1/12 

1/15 


1/60 



(ID 


( 12 ) 
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will give the equations ( 6* ) through ( 16* ) (except ( 14* ) 5 ) in Table I in the 
following order: 6 


(6*) 

(7*) 

(8*) 

(9*) 

(10*) 

- 

(11*) 

- 

- 

(12*) 

- 

- 

(13*) 

(15*) 

- 

( 14a* ) # 

- 

- 

(14b*)# 

- 

- 

(16*) 

- 

- 


These 7 together with equations (5) and (6) are the Runge-Kutta 
equations for the fifth order with N = 5. (See footnote 1. ) 

REDUCTION OF THE RUNGE-KUTTA EQUATIONS BY SUBSTITUTIONS 


The solutions of the Runge-Kutta equations are effected by arbitrarily 
assuming (except for some restrictions which will be noted later) the values 
of the parameters a . and with these values solving the first five equations in 

Table I for (or, in the more general case, the (m + 1) equations given by 

( 5) ) . Since this procedure maybe universally used for all orders and for either 
exact or approximate solutions, the special technique which is applicable to 
this system of equation is outlined in Appendix A . 


Equation ( 14* ) is obtained by adding ( 14a* ) and ( 14b* ) (the equations 
marked # in ( 13) ) . 

®It will be observed that the numbers of the equations in Table I are 
distinguished by astrisks. 

Equations which are not relevant to our process are indicated by dashes 
in the equation list ( 13) . 


8 



With a. and co^ known, the remaining Runge-Kutta equations (the equa- 
tions ( 6) and the equations from the appropriate matrix equation analogous to 
equation ( 12) ) are then used to solve for /3„ . It is, therefore, to be noted 

that any transformation which converts the latter Runge-Kutta equations into 
those given by equation ( 5) is especially valuable, since this reduces the 
number of equations in the system which determine the parameters. 

We will now show how, by some simple substitutions, we can substanti- 
ally reduce the number of Runge-Kutta equations in Table I. 

In the first column of the matrix ( 11) , let us set 


(i) _ 

c 3 

@32 a 2 

= <*1/2 

(14a) 

<0 

C 4 

P42&2 + ^ 43^3 

<N 

3 

II 

(14b) 

0) 

c 5 

042«2 + ^ 53«3 + PtoPi 

= a !/ 2 = 1/2 

(14c) 


It is immediately evident that equations ( 6* ) , ( 9* ) , ( 11* ) , and ( 12* ) 
in ( 13 ) will be transformed respectively into ( 3 * ) , ( 4 * ) , ( 5* ) and ( 5* ) by 
these substitutions, provided w 2 = 0, so that the former equations may be 
removed from further consideration. 8 

The display in the equation list ( 13) will now leave effectively the 
following equations: 


- 

( 7 *) 

(8*) 

- 

(10*) 

- 

(13*) 

(15*) 

- 

( 14a* ) 

- 

- 

( 14b* ) 

- 

- 

(16*) 

- 

- 


(15) 


8 For higher orders we can also make similar substitutions for the 

columns in (11); e. g. , we can let cj 2 ^ = c^/3 (i = 3, 4, 5 ) where c^ is 

defined by equations ( 9) . 
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Using the notation in matrices (7) and ( 8) , we can write the following 
equations on the assumption that equation ( 14) is used in the matrix C ( 11) : 


<°) 2 . (0) 2 . (0) 2 
72 “2 + 73 « 3 + 74 <4 

= 1/12 

(7*) 

( 0) 3 , (0) 3 , (0) 3 

72 a 2 + 73 <4 + 74 «4 

= 1/20 

(8*) 

(l) 2 1 (l) 2 j- ( i) 2 

72 “2 + 73 «3 + 74 4 

= 1/15 

(10*) 

(0) 2 . (°) 2 
73 4 + 74 4 

= 1/12 

(13*) 

(0) 3 , (0) 3 

73 “3 + 74 4 

= 1/20 

( 14a* ) 

(l) 2 a. (l) 2 
73 4 + 74 4 

= 1/15 

( 14b* ) 


Let us now assume, for the fifth-order approximation, that 72 = 0 

and furthermore, that 72 ^ Is approximately zero; or, sufficiently small so 

that we may neglect it. 9 It is immediately apparent from equations ( 16) that 
each equation can be paired and that another reduction in the number of 
equations that have to be satisfied has been achieved. Thus equation list (15) 
may now be replaced by 


- 

(13*) 

( i4a* ) 

- 

(14b*) 

- 

(13*) 

(15*) 

- 

( 14a* ) 

- 

- 

(14b*) 

- 

- 

(16*) 

- 

- 


Our problem is to determine the parameters /3„ from these equations. 

Something about the nature of the necessary approximations to solve these 
equations may be obtained from the following analysis. 


3 It should be observed that the argument presented here should hold for 
higher orders as the equations ( 16) will be unchanged except that there will 
be more terms. 
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It may readily be verified by inspection that the equations ( 13* ) , 
( 14* ) , 10 ( 15* ) and ( 16* ) may be represented by the matrix equation 


/ 1 1 1 ‘\ 

/ n>4o: 2 (^32/^43^ 


/l/24\ 

[(0:3+0! 4) ( 0:3+0: 5 ) («4+0! 5 ) ( 0! 4+0:5) 

^5^2(^32^53) 


f ' 

7/20 

i Q!2 a 2 0 i 2 <x z 

^ 5^2 (^42^54) 


1 1/60 

Vo 0 0 012/352 J 

\^ 5 of 3 (^ 43/ 3 54 )/ 


\ 1/120/ 


The (4x4) matrix in (18) is singular so that the determinant of the 
homogenous system is zero. Therefore, in order that the system (18) have 
solutions. Hie augmented matrix for the non-homogenous system must be of 
rank three, in which case it may be verified that the following relation must 
hold 


( g -3 - Q ! 2 ) a z + P32&2 ( 5 «2 -2) = 0 


(19) 


provided the roots a z , a 4, and P's are distinct. 

Substituting equation (14a) into equation (19) , we get 

0:2(5 o 3 -2) = 0, (20) 

which together with tiie condition arising from the choice u> 2 = 0 (see ( 9- A) , 
Appendix A) restricts the choice of o. when o 5 = 1 and u> 2 = 0. 

Neither of the two solutions that satisfy equation (20) is tenable. o 2 can- 
not be zero, nor o 3 = 2/5, since by the condition (11) in Appendix A it is associ- 
ated with 04 = 1 ; and we should recall that a repetitive o is excluded. 

We cannot, therefore, satisfy equation (20) exactly; but we can get an 
approximate solution if we make o: 2 small. We shall proceed on this supposition. 

10 See footnote 5, page 8. 
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OTHER TRANSFORMATIONS 


We may now impose certain additional conditions on our parameters 
and obtain further significant reductions in the Runge-Kutta equations. Hufa 
showed [6, p. 204; 5, p. 8] that what is essentially in our notation, 


u. 

y 1 , = ^ • ( ^ — » j~ N ; P — 0 > j — i , 

j _ j 1 X J i J X J 


will hold if 



u >2 = 0 and 



( 21 ) 


and, therefore, as may be seen by using matrix (8) 11 

y. (0) = w (1 - a.) , j = 3,4 N . (22) 

With these relations it is immediately evident with the help of equa- 
tions (16) that the equations (7*), (8*), (13*) and (14a*) are satisfied; i. e. , 
these equations are converted into differences of equations in equation ( 5) . 
However, the importance of these transformations lies not only in these simpli- 
fications but also in the fact that they give us simple linear relations from which 
to calculate the in place of the non-linear Runge-Kutta equations. These 
linear relations are obtained simply by equating the two definitions of 

y! the one given by matrix ( 8) and the other by equation ( 22) . 


A special , but important, example of this is the relation 
y N-l = ^N^N.N-l 


(23) 


which follows easily from the definition in matrix ( 8) ; but it must be noted it 


n See Appendix D for establishing the validity of these transformations. 
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holds only in this simple form for the ( N-N) case. Combining now the two 
evaluations for given by equations (22) and (23) , we get 


Y N-1 W N^N,N-i W N-i (1 “n-1* 


( 24 ) 


This is an important relation as it immediately gives us one of the parameters 
^N, N-i* 


THE PARAMETERS FOR THE N™ ORDER WITH N EVALUATIONS 


It might be well now to examine the number of equations that are 
available for finding approximate Runge-Kutta solutions for some of the higher 



process can readily be seen from equations (2) to be equal to (N-i) (N-2)/2. 
The number of equations available for finding these parameters are shown in 
Table II. When the number of available equations ( given in the last column of 
the table) is less than the number of parameters /3„ to be found, we can make 

up this deficiency by using some of the Runge-Kutta equations. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administrations 
Huntsville, Alabama, April 10, 1967 
039-00-24-00-00 


12 When j = 1 , 0 
have been calculated. 


can be found by using equation (6) 


after the other /?'s 
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TABLE II. THE NUMBER OF EQUATIONS AVAILABLE FOR FINDING THE PARAMETERS 0 (j # 1) 


Order , 

No of Parameters 
/3 (j ^ 1) Required 

= (N-i) (N-2)/2 

cj‘Wj/2 
is (N-2) 

Number of equations available from: 

• « = 0!®/3 
is (N-2) * 

(21) and (22) 
is (N-2) 

sU 

»-*- 

tl 

o 

Total 

Available 

Eqs. 

5 

6 

3 

0 

3 

0 ** 

6 

6 

10 

4 

0 

4 

i 

9 

7 

15 

5 

0 

5 

i 

11 

8 

21 

6 

6 

6 

i 

19 

9 

28 

7 

7 

7 

i 

22 


* 


When used 

72 is assumed to be approximately zero in this case. 




APPENDIX A 

THE WEIGHTS Uj IN TERMS OF o j 


The weights co. in equation ( 5) can be found in terms of ce. by an effec- 
tively simple method [11]. Though the method is general we shall demonstrate 
this technique for the fifth-order cases with five evaluations. 1 

Construct the polynomial 


F(ai) = (a - aj) (a - a 2 ) (a - a 5 ) 

= a 5 + ajG! 4 + a 2 Q! 3 + a 3 a 2 + a 4 a< + a 5 


(1-A) 


where the ct^ are parameters in our system of equations (a l = 0) . 


We can then construct a second polynomial 


G(g 0 = g 4 a 4 + g 3 a! 3 + g 2 Q! 2 + gi a + g 0 


(2-A) 


whose coefficients are easily determined from the coefficients in equation (1) 
and the right hand members in the matrix (5) . Then, the weights in the system 
of equations in matrix ( 5) are given by [11]. 


G(a.) 
W i = F'(a.) 


i = 1,2, 3, 4, 5. 


(3-A) 


We shall limit ourselves for the present to the particular set of solutions 
where a 5 = 1 and 002 = 0, so that all subsequent formulas are specialized to these 
values. We should also note that, although cxi does not explieity occur in the 
Runge-Kutta equations , since it is equal to zero , it does appear in polynomials 
(1-A) and (3-A) . 


4 This is the case where N = m + 1 in matrix (5) . If , as is usually the 

case, N > m+ 1 we can revert to the foregoing case by selecting N - (m + 1) 

values of o>.. 

x 
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Hence, when a 1 - 0 and a 3 = 1, ( i-A) takes the special form 

F(a) = a 5 + a t a 4 + a^a 3 + a 3 a! 2 - (l + a 1 + a 2 + a 3 )a (4-A) 

where , since the coefficients in ( 3-A) are the symmetric functions involving its 
roots, 


a t = ~(a 2 + a 3 + q ! 4 + 1 ) 

a 2 = (oi 2 oi 3 + a 2 a i + + (ot 2 + a 3 + a 4 ) (5-A) 

a 3 = -a 2 a 3 a 4 - ( a 2 a 3 + a 2 a 4 + a 3 a! 4 ) 


while the coefficients of G (a) in (2-A) are given by [ii] 


Si = i 

g 3 = a 4 + i/2 

g 2 = ^ + (i/2)a 4 +1/3 (6-A) 

gi = a 3 + ( i/2) a 2 + (1/ 3)a 4 +1/4 

So ~ ( - l/2) a 3 ~ ( 2/3) a^ - ( 3/ 4) aj - 4/5 

We shall also find useful the cubic derived from equation ( 4-A) 
f(o!) = a 3 + (a 4 +1) o; 2 + (a 2 + a 4 + 1) a + (a 3 + a^ + a 4 + 1) = 0 (7-A) 

which has the roots a 2 , a 3 and Q! 4 . 
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The Condition For ^ = 0 


Since we are assuming that o >2 = 0, an important condition may be 
derived when we note from equation (3-A) that G ( a 2 ) = 0 must hold, provided 
F* (a 2 ) =£ 0. However, we can exclude F' (a 2 ) = 0, since this would denote a 
double root a 2 in F (a). We can expedite the determination of this condition if 
we evaluate instead 


2G(a 2 ) - (20-2 + 1) f (« 2 ) = 0 


( 8-A) 


since by equation (7-A) 0:2 is also a root of f (a) =0. 

After somewhat tedious algebraic operations on substituting the a's 
given in equation ( 5-A) into equation ( 8-A) , we get 


100:3 a i - 5(03 + 0:4) +3 = 0. 


(9-A) 


This is a relation that must subsist between the two parameters when 
o >2 = 0, ffj = 0 and a 5 - i. 

Formulae For The Weights Uj 


Since G ( qj 2 ) = 0, we can write equation (3-A) in the form 
(a - a 2 ) h (a ) 

M i = ~F»Ta ^) ’ i = 1. 2, 3,4, 5 (10-A) 

where, by dividing G {a) synthetically by (a - a 2 ) , or factoring by other means, 
we find 


h( cv) = a 3 + hia 2 + h 2 a +h 3 (ii-A) 
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where 

hi = « 2 + g 3 = - (“3 + Ci 4 + 1/2) 

h 2 = h^j + g 2 = a 3 Q! 4 + ( 1/2) (<* 3 + «4) - i/6 

( 12- A) 

h 3 = h 2 o; 2 + gi = - a 2 ot z a 4 - (1/2) (<x 2 a z + + a 3 o; 4 ) 

+ (1/6) (a 2 + cx 3 + 0-4) - 1/12 . 

Using equation ( 9-A) , these may be simplified so that we get 

h(oi) = a 3 - (q ! 3 + a 4 + i/2)o! 2 + (a 3 + a 4 - 7/l5)a ( 13-A) 

- (1/60) [5(a 3 + a 4 ) - 4] . 

It will be observed that the coefficients of h (a) in equation ( 13-A) are 
independent of a 2 . 

With the aid of the first equation in ( i-A) we can readily find the 
F* (a.). 

These are: 

F* (a-i) = ch 2 q; 3 q ;4 , = 0 

F’ (a 2 ) = a 2 (a 2 - a 3 ) (a 2 - 0-4) ( a 2 - 1) , 

F* (o; 3 ) = a 3 (a 3 - a 2 ) (a 3 - £*4) (a 3 - 1) , 

( 14- A) 

F* (a 4 ) = a i (a i - a 2 ) (0:4 - a 3 ) (0:4 - 1) , 

F 1 (ce 5 ) = Qi5(a! 5 - a 2 ) (a 5 - a 3 ) (a 5 - a 4 ) 

= ( 1 - a 2 ) ( 1 - a 3 ) ( 1 - a 4 ) , Q! 5 = 1 . 

Using these in ( 10-A) , the formulae for the weights, we get 

co 1 = — h ( 0) /ol 3 ol^ (15— A) 
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o 

II 

CM 

3 

( 16-A) 

co 3 = Mag) /a 3 (oi 3 - Q! 4 ) (ag - 1) 

(17-A) 

= h(Q! 4 )/Q!4(Q! 4 - CKg) (Q! 4 - i) 

( 18- A) 

co 5 = h(Q! 5 )/(Q!3 - 1) (a 4 -i) , a 5 = 1 

(19-A) 


It should be observed that, in view of our previous observation on the 
coefficients of h (a) in equation (13-A) , the also are independent of oi 2 . 

An Alternate Method 


An alternate method (which we will find useful to refer to later) is to 
solve the linear system of equations (5) by determinants. We will illustrate this 
for the sixth-order case (N = 6, m = 5) as our method may easily be extended to 
higher orders. 

For u >2 we have the ratio of the two determinants: 


I 


0 

1/2 

a 3 

“4 

“5 

a 

0 

1/3 

a 3 

a\ 

a 5 

a 

0 

1/4 

« 3 

« 4 

a\ 

a 

0 

1/5 

a 3 


a 5 

a 

0 

1/6 

« 1 

“4 

«5 5 

a 


(20-A) 
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and 


0 a 2 a 3 

0 «2 a 3 

0 al “3 

0 a| al 

0 al a® 


Oi 4 

«5 

“6 

a\ 



a\ 

“5 

“6 

a\ 

« 5 

o\ 

Ci\ 

a 5 

<4 


(2i-A) 


These determinants may be further simplified for we can reduce the 
order of determinants (20-A) and (21-A) by taking the cofactor of 1 in the first 
column. Then if we set a 6 = 1 and subtract adjacent rows, we get for determinant 
(20-A) 

1/6 a 3 ( t-a 3 ) 0 ' 4 (l- 0 ' 4 ) a 5 (l-o; 5 ) 0 

1/12 a\{ i~o? 3 ) a\( i-a 4 ) a\{i-a 5 ) 0 

1/20 a^U-ag) o;|(i-a 4 ) Q! 5 3 (l-Q! 5 ) 0 (22-A) 

1/30 a!|(l-a3) q;|( 1-Q! 4 ) a|(l-a 5 ) 0 

1/6 af a\ 

which reduces to 

1/16 a 3 a 4 

1/12 a 3 a\ 

(i-a 3 ) (l-o; 4 ) (l-a 5 ) 

1/20 a\ a\ 

1/30 a 3 a 4 
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Similarly determinant (21-A) can be reduced to 


(1— G! 2 ) (l-a 3 ) ( 1-0(4) (l-«s) 


«2 

CO 

3 

« 4 

a 

«2 

“3 

a\ 

a' 

COM 

3 

“3 

a\ 

a' 

a 2 


a i 

a 


( 24-A) 


so that o >2 ( l-« 2 ) is given by the ratio of the determinants in (23-A) and 
(24-A): 


1/16 

a 3 

a 4 

«5 


«2 

«3 

«4 

“5 

1/12 

a 3 

a\ 

«5 


«2 

“3 


a 5 

1/20 

a 3 

“4 




«3 

«4 

“5 

1/30 

4 

« 3 

«4 

a i 


4 

“2 

“3 

“4 

“5 


Similar ratios may be found for the expressions o^(i - a ) . 

It is evident that the vanishing of the determinant (23-A) is a condition 
for co 2 = 0 provided a 6 = 1 and a # 1 (i = 3, 4, 5 ). 2 


2 To obtain this condition for the fifth-order case delete the last row and 
column in the determinant (21-A) and set Q ! 5 = i. The vanishing of this deter- 
minant is the relation ( 9) . 
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APPENDIX B 

THE PARAMETERS FOR THE FIFTH ORDER WITH FIVE EVALUATIONS 


We shall now show how the six parameters (j + 1 ) 1 for the approxi- 
mate fifth-order Runge-Kutta process with five evaluations may be found. We 
assume that the values of a 3 and a 4 have been selected (subjected to the condi- 
tion ( 9-A) ) and that the u>'s have been obtained as indicated in Appendix A. 

We can now proceed as follows — expressing the parameters in terms of the 
known a 3 , 0-4 and o>. : 2 


i) Calculate /3 32 by equation ( 14a) : 
032 = <2 3 / 20*2 


2) Calculate jS 54 by using equation (24) which defines y\ in two ways 



- w 5054 


^(i-C^) , 


so that 


054 = a> 4 ( 1-0:4) /( t > 5 


3) Calculate /? 43 using equation ( 16* ) : 


043 “ 


60 CO 50:3 £54 


60 Q! 3 C04( 1-0-4) 


*When j = 1, /9„ can be found by using equation ( 6 ) after the other (3's 
have been calculated. 


z As previously observed the ccl are independent of a 2 so that in the sub- 
sequent results a 2 is a free parameter. 
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4) Calculate 04J. using equation ( 14b) : 


042 = 




( 0 ) 

5) Calculate 0 52 using y 2 = 0: 


052 — 


_1_ 

^5 


(^ 303 2 + 0 ) 4043 ) • 


6) Calculate 0 53 using equation ( 14c) : 
053 Q! 3 = 2 ” ^52^2 + 054^4 ) 

or, instead use 

73 ^ = 0)4043 + ^5043 = w 3 ^ ^ -a 3 ) 


the two results should check. 



APPENDIX C 

THE SIXTH-ORDER RUNGE-KUTTA APPROXIMATION 
WITH SIX EVALUATIONS 


We shall outline a step-by-step procedure for obtaining the parameters 
for the sixth-order Runge-Kutta approximation with six evaluations. 

As outlined in Appendix B for the fifth-order case, we can in a similar 
way select values for (with some restrictions) and find the u. . As we 

assume equations (21) to hold, it is true that otg = 1 and u> 2 = 0; and it is this 
last relation which dictates a condition between the o:. (i = 3, 4, 5) analogous 

to ( 9- A) , and which we will not stop to derive. 

As we may see from Table II we can assume that the following hold: the 
four relations analogous to equation (14) ;* the four relations in equations (21) 

and (22) ; and finally, that y2 ^ = 0; 2 a total of nine relations for ten unknown 
parameters. We are thus free to make use of another Runge-Kutta equation; 
and we might make use of the equation analogous to equation ( 16* ) which was 
found useful in the fifth-order case. This is: 3 

w 6^32 ^43^54^65 a 2 = ‘ (i“C) 


We now proceed as follows: 

1) Calculate fi 32 using equation ( 14a) : 


^32 — «l/2«2 • (2-C) 


1 C i (1) = a i^ 2 (i = 3 ’ 4 ’ 5 > 6 >- 

2 See footnote 9 . Since the equations in ( 16) will hold for the sixth-order 
case also (except for additional terms) , it may be seen that in this group ( 10*) 

would revert to equation ( 14b* ) if we set y 2 ^ = 0. This assumption is vital to 
the simplifications necessary to reduce the complexity of the computations 
given here. 

3 See footnote 1, page 3. 
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2) Using the relations yg = co e /S 65 = co 5 ( l-a 5 ) given by equation (24) : 


^65 = co 5 (i-a! 5 )/a> 6 


(3-C) 


(The subsequent /3’s are all calculated in terms of the unknown /3 54 .) 
3) Using equation ( i-C) calculate : 


^ (. 6 J ujrPrr(x\ ) Pra 


(4-C) 


4) Calculate /3 42 using equation (14b) 


~ a 2 [if " ^ 43 “ 3 1 


(5-C) 


5) Solve for /3 52 using the relations: 


(o) 

72 ~ Uz&7& + Ci) 4 ^ 42 + Wg£g 2 + COg/362 - 0 

72 ^ = (^Z^22 a Z + U4P42&4 + ^ 5 ^ 52^5 + w 6 ^ 62 a 6 = 0 


( 6-C) 


which by subtraction gives (Since a 6 = i) 


72 ^“72 ^ - w 3< ^■~ 0i z)^Z2 + w 4<^-~ Q! 4)^42 + w 5< i~ a $)Pz2 ~ 0 • (7-C) 


6) Solve (in terms of /3 54 ) for j3 53 using 


C 5^ = / 5 52«2 +PS3&3 + Pu a 4 = a \/ 2 


( 8-C) 
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7) Substitute the foregoing values of /3. . into equations (21) and (22) 

where > 73 and y \ are expanded by equatibns ( 8 ) and find /S^ and 
/ 8 g 4 in terms of 


8 ) Finally determine j3 64 by substituting the above values of /3„ into 



- ^ 62 Q! 2 + 063 a 3 + + ^ 65 a 5 ~ a \/% ~ 1/2 


(9-C) 
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APPENDIX D 

THE VALIDITY OF THE TRANSFORMATION 5; * (I - a jl 


It would be desirable to examine the validity of assuming equations (21) 
and (22) when seeking approximate Runge-Kutta solutions. We will illustrate 
this with an analysis for the sixth-order case with six evaluations as the pro- 
cedure can easily be adapted to higher order cases. 

Let us examine the effect of these transformations on the linear 
system analogous to (6*) , (7*) and (8*) — the equations in the first row in 
the equation list ( 13) . (See also equations ( 16) . ) For the sixth-order there 
will be an additional equation in the system and each equation will have an addi- 
tional term; these will be given by (see footnote 1 page 3 , Table I) 



n 

a 2 



n 

“3 



n 

«4 



n 
a 5 


= l/(n +1) (n + 2) 


n = 1, 2, 3, 4 . 


(1-D) 


It is immediately evident that the substitutions of equations (21) and (22) 
into the system (1-D) gives us simple combinations of equations in system (5) 
if a 6 = 1, so that the advantage of these substitutions is apparent. However, we 
must see if equations (21) and (22) are acceptable as solutions of the system 
(1-D). 

If we solve the linear system (1-D) for yjj using determinants, we get 


1/6 

“3 

“4 

a 5 


“2 

“3 

«4 

“5 

1/12 

a l 

“ 4 

“5 





«5 

1/20 

“3 


“5 

• 

al 

« 3 

a l 

a l 

1/30 

«: 3 

a i 

<4 


<4 

<4 

a i 

a 5 


which we have shown in equation (25-A) to be a> 2 ( l-a^) • 

In the same way we can quite readily justify the other transformations 
in equation ( 22) . 
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